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Interface dynamics of a metastable 
mass-conserving spatially extended diffusion 
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Abstract 

We study the metastable dynamics of a discretised version of the mass-conserving 
stochastic Allen-Cahn equation. Consider a periodic one-dimensional lattice with N 
sites, and attach to each site a real-valued variable, which can be interpreted as a spin, as 
the concentration of one type of metal in an alloy, or as a particle density. Each of these 
variables is subjected to a local force deriving from a symmetric double-well potential, 
to a weak ferromagnetic coupling with its nearest neighbours, and to independent white 
noise. In addition, the dynamics is constrained to have constant total magnetisation 
or mass. Using tools from the theory of metastable diffusion processes, we show that 
the long-term dynamics of this system is similar to a Kawasaki-type exchange dynam¬ 
ics, and determine explicit expressions for its transition probabilities. This allows us 
to describe the system in terms of the dynamics of its interfaces, and to compute an 
Eyring-Kramers formula for its spectral gap. In particular, we obtain that the spectral 
gap scales like the inverse system size squared. 
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1 Introduction 

The low-temperature dynamics of spatially extended systems often displays metastability: 
these systems can spend considerable amounts of time in configurations that have higher 
energy than their ground state. Well-known examples of such phenomena are supercooled 
water, which remains liquid at temperatures below 0°C, a supersaturated gas, which does 
not condensate although this would be thermodynamically more favourable, and a wrongly 
magnetised ferromagnet. 

Much research effort has been dedicated to the study of metastable lattice systems, 
such as the Ising model at low temperature. This has led to very precise results on the 
time the system spends in metastable equilibrium, on the way it moves from a metastable 
to a stable state by creating a critical droplet, and on the shape of this droplet. See for 
instance m for a review on Ising models with Glauber (spin flip) dynamics and lattice gases 
with Kawasaki (particle/hole exchange) dynamics, and [27] for results based on the theory 
of large deviations. A considerably more difficult case arises when there is no underlying 
lattice given a priori, but particles instead evolve in M d , and one wants to describe processes 
such as crystallisation. For recent results in this direction, see for instance mmm- 

Another type of models whose metastable behaviour is understood in detail are diffusion 
processes described by stochastic differential equations with weak noise. A general large- 
deviation approach to these equations goes back to the work of Freidlin and Wentzell [20], 
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which provides many results on transition times between attractors and on the long-time 
dynamics. In the case of reversible diffusions (that is, those satisfying a detailed balance 
condition), metastable timescales are governed by the so-called Eyring-Kramers formula, 
derived heuristically in |18l 25] , and first proved in a mathematically rigorous way in maim. 
See for instance [3] for a recent survey on various methods of proof and extensions of the 
result. 

A spatially extended system of coupled diffusions, which can be considered of interme¬ 
diate difficulty between lattice systems with discrete spins and systems of particles evolving 
in M rf , was introduced in mm- In this model, the spins are still attached to a lattice 
(which is periodic and one-dimensional of size N), but they take values in M instead of 
{—1, +1}. Each spin feels a local symmetric double-well potential with minima in ±1, and 
is coupled ferromagnetically to its nearest neighbours. In addition, each spin is subjected 
to independent white noise. For weak coupling, the dynamics of this system was shown to 
be similar to that of an Ising model with Glauber spin-flip dynamics. Indeed, the energy 
of configurations increases with the number of interfaces , defined as pairs of neighbouring 
spins having different sign. As a consequence, the system favours configurations with few 
clusters of spins having the same sign. On the other hand, when the coupling scales like 
N 2 , the system converges as N —> oo to an Allen-Cahn SPDE with space-time white noise, 
whose metastable behaviour was studied in m eg. 

A natural question that arises is whether one can construct a similar system, with con¬ 
tinuous spins attached to a discrete lattice, but whose dynamics for weak coupling resembles 
Kawasaki exchange dynamics instead of Glauber spin-flip dynamics. In other words, one 
would like to impose that the total magnetisation (or the total mass in lattice gas termi¬ 
nology) is conserved. A simple way of doing this is to start with the potential energy of 
the system considered in and to constrain it to the hypersurface where the sum of 

all spins is constant, say equal to zero. This is nothing but the discretised version of the 
mass-conserving Allen-Cahn equation introduced in |29| . The objective of the present work 
is to study the metastable dynamics of this model. 

It is quite easy to see that in the uncoupled limit, the potential energy of the constrained 
system is minimal when exactly half the sites have value +1, while the other half have 
value —1. Such states have a clear particle system interpretation: just consider each +1 as 
a particle and each —1 as a hole. As in the unconstrained case, for weak positive coupling, 
the energy of configurations increases with the number of interfaces. Therefore the ground 
state consists of the configurations having exactly one cluster of particles and one cluster 
of holes, separated by two interfaces. Higher-energy configurations have more clusters and 
more interfaces. Thus if the system starts in an excited state with many interfaces, one 
expects that its clusters will gradually merge, reducing the number of interfaces, until the 
ground state is reached (Figure [I]). 

While our analysis will show that this picture is essentially correct, there is a complica¬ 
tion due to the fact that particle/hole configurations are not the only local minima of the 
potential energy. Somewhat unexpectedly, there turn out to be many more “spurious” local 
minima, whose coordinates are not close to ±1. The way around this difficulty is to realise 
that all spurious configurations have a higher energy than the particle/hole configurations. 
Therefore the long-term dynamics will spend most of the time near the particle/hole config¬ 
urations, with occasional transitions between them. Our main result is the characterisation 
of this effective dynamics. 

This paper is organised as follows. In Section [2j we give a precise definition of the 
considered model. In Section [3j we describe the potential landscape of the model, meaning 
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Figure 1. Example of evolution of the constrained system (2.7) with N = 512 sites. Space 
goes from left to right, and time from top to bottom. Blue and red correspond to spin values 
close to —1 and 1 respectively. The system starts in a configuration with 40 interfaces, many 
of which disappear quickly. At the end of the simulation, the number of interfaces has been 
reduced to 4. Parameter values are £ = 0.02 and 7 = 16. This coupling intensity, which 
is much larger than considered in this work, has been chosen to obtain transitions on an 
observable timescale. 


that we find all local minima of the potential energy, and describe how they are connected 
by saddles with one unstable direction. Section [4] uses the notion of metastable hierarchy to 
show that the dynamics indeed concentrates on particle/hole configurations, and derives the 
effective dynamics on these states. In Section [5] we use this information to characterise the 
evolution of interfaces, and we derive a sharp estimate for the spectral gap of the system, 
which determines the relaxation time to equilibrium. Section [ 6 ] contains concluding remarks, 
while most proofs are postponed to the appendix. 

Notations: If i ^ j are integers, \i,j\ denotes the set {i,i + 1,... ,j}. The cardinality of 
a finite set A is denoted by |A|, and A = B C indicates that A = B U C with B and C 
disjoint. We write 1a for the indicator function of the set A , l n or simply 1 for the identity 
matrix of size n x n, and 1 for a column vector with all components equal to 1. Finally, 
we write E^[-] for expectations with respect to the law of the diffusion process started with 
distribution /q and E x [-] in case q is concentrated in a single point x. 

Acknowledgements: The idea of studying the constrained process considered in this work 
goes back to a question Erwin Bolthausen asked after a talk given in Zurich by the first 
author on the unconstrained model studied in mm- 

2 Definition of the model 

Consider the potential —> M defined by 

N N 

V ' y ( x ) = J 2 U ( Xi ') + - 4 : J 2( Xi + 1 ~ Xi ') 2 ’ ^(0 = ^ 4 -^ 2 , ( 2 - 1 ) 

2=1 2=1 

where N ^ 2 is an integer and 7 ^ 0 is a coupling parameter. We also make the identification 
xyr+i = xi, that is, we consider periodic boundary conditions. Thus x can be considered 
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either as an element of W N , or as an element of M A , where A is the periodic lattice Z/JVZ. 

The potential allows to define a diffusion process by the stochastic differential equa¬ 
tion 

d xt = —’VVy(xt) d t + V2edWt , (2-2) 

where Wt is an A r -dimensional Wiener process, and £ ^ 0 is a small parameter measuring 
noise intensity. The dynamics of this system has already been studied in PE]- Here we are 
interested in a different system, obtained by constraining the diffusion to the hyperplane 

^Xi = o| . (2.3) 

^ »=l J 

To define its dynamics, let R be an orthogonal matrix mapping the unit normal vector to S 
to the Nth. canonical basis vector e^- Let Vy(y) = V y (R~ 1 y), and define the dynamics by 

d yi,t =- 7 —- df + V2e dW ht , i = 1,..., N - 1 , 

oyi,t 

VN,t = 0 , (2.4) 


where Wij, ■ • •, ITv-ir are independent Brownian motions. Then xt is by definition the 
process xt = R~ 1 yt■ It is easy to check that this definition does not depend on the choice 
of R. 

An equivalent way of defining the dynamics is to write 

d Vt = [-VH 7 (y t ) + (VVy(y t ),e N )e N \ + V2edW t (2.5) 


where Wt is an (N — l)-dimensional Wiener process. Indeed, the extra term precisely 
ensures that the A r -th component of the drift term vanishes. Transforming back, we obtain 
the equation 


d x t = 


-w 7 (x t ) + -(vy 7 (x t ),i)i 


df + \/ 2 i d Wt 


t •) 


( 2 . 6 ) 


where 1 denotes the vector with all components equal to 1 (hence the normalisation 1/N), 
and Wt = R~ 1 Wt is a Brownian motion on S. When written in components, the resulting 
dynamics takes the form 


dx i)t = 


i N 

j- % 

+ - (Xi+i,t - 2 Xij + x»_i jt ) - — 2_^ 


j = 1 


df + V2e d W, 




(2.7) 


where /(^) = —U'(£) = £ — (and the Wj.t are no longer independent). Note that this is 
a discretised version of the mass-conserving Allen-Cahn SPDE 

d t u(t,x) = ~/Au(t,x) + f(u(t,x)) - f(u(t,y))dy + V2s€(t,x) (2.8) 


with space-time white noise £ on S. The nonlocal integral term indeed ensures that the total 
mass u(t, x ) dx is conserved. This equation was introduced in |29j in the case without 
noise, and considered recently in [ 1 ] in the case with noise. 


Systems of the form ( 2 . 2 ) admit a unique invariant probability measure with density 


y(x) = ^e- v ^, 


(2.9) 
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where Z is the normalisation constant, and are reversible with respect to /r. Analogous 
statements hold true for the system constructed here (except that [i is concentrated on the 
hyperplane S). The questions we thus ask are the following: 

• How long does the system take to relax to equilibrium? 

• What are the typical paths taken to achieve equilibrium, when starting in an atypical 
configuration? 

• Can the system be approximated by a coarse-grained process visiting only local minima 
of the potential? What does this coarse-grained process look like? 


3 Potential landscape 

3.1 The transition graph 

For a general system of the form ( |2.2[ ), let 

5 = {x e R n : W 7 (x) = 0} (3.1) 


be the set of all stationary points of Fy. A stationary point x* £ S is called non-degenerate 
if its Hessian matrix V 2 Fy(x*) has a nonzero determinant. We will assume for simplicity 
that all stationary points of Fy are nondegenerate (see however [BJ for results on systems 
with degenerate stationary points). 

The Morse index of a nondegenerate stationary point x* is the number of negative 
eigenvalues of the Hessian V 2 Fy(x*) (i.e., the number of directions in which Fy decreases 
near x*). For each k £ [[0,1V]], let S} ; denote the set of stationary points of index k. The 
set So of local minima of Fy and the set 5i of saddles of index 1 (or 1 -saddles) are the most 
important for the stochastic dynamics for small e. 

By the stable manifold theorem, each 1-saddle has a one-dimensional unstable manifold 
consisting in two connected components. Along each component, the value of Fy has to 
decrease, and therefore (since Fy is confining) both components have to converge to a local 
minimum of Fy. Let Q = (So,£) be the unoriented graph in which two elements of So are 
connected by an edge in £ if and only if there exists a 1 -saddle z £ <Si whose unstable 
manifold converges to these local minima. 

Roughly speaking, the stochastic system behaves for small noise intensity e like a Marko¬ 
vian jump process (or continuous-time Markov chain) on So, with jump rates related to the 
potential differences between local minima and 1-saddles. This is the basic idea imple¬ 
mented in m Chapter 6 ], and there are many refinements on which we will comment in 
more detail below. 

In the case of the potential (2.1) without constraint, the potential landscape has been 
analysed in [ 6 j. In particular, the following properties have been obtained: 

• If 7 = 0, the set of stationary points is given by S = {—1,0,1}^. The local minima 
are given by So = {— 1 , 1 }^ and the 1 -saddles are those stationary points that have 
exactly one coordinate equal to 0. They connect the local minima obtained by replacing 
the 0 coordinate by —1 or +1. Thus the graph Q is an IV-dimensional hypercube, 
with transitions consisting in the reversal of the sign of one coordinate, which can be 
interpreted as spin flips. 

• There exists a critical coupling 7 * (IV), satisfying 7 * (IV) ^ j for all N, such that the 
transition graph Q is the same for all 7 £ ] 0 , 7 *(IV)). Thus the local minima and 
allowed transitions are the same for weak positive coupling as in the uncoupled case. 
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Figure 2. Transition graph of the unconstrained system for N = 4 and 7 = 0. Black and 
white circles represent respectively coordinates equal to 1 and to —1. The two configurations 
( 1 ,- 1 , 1 ,- 1 ) and (— 1 , 1 ,— 1 , 1 ) are not shown, because they correspond to non-optimal 
transitions as soon as 7 > 0 . 


What changes, however, is that some transitions are easier than others when 7 > 0: the 
systems prefers transitions that minimise the number of interfaces, that is, the number 
of nearest neighbours with a different sign (Figure [2|. The stochastic dynamics is thus 
very close to the one of an Ising model with Glauber spin-flip dynamics. 

• For 7 increasing beyond 7 * (IV), the system undergoes a number of bifurcations that 
reduce the number of stationary points. In particular, for 7 > 1/(2 sin 2 ( 7 r/iV)) the 
system synchronises: there are only two local minima given by ±( 1 , 1 ,..., 1 ), connected 
by the only 1 -saddle which is at the origin. 

Our aim is now to obtain similar results for the graph Q of the constrained system, 
starting with the uncoupled case 7 = 0 and then moving to small positive 7 . 


3.2 The uncoupled case 


We consider in this section the dynamics of the constrained system in the uncoupled case 
7 = 0 . The above definitions of Sq, Si and Q can be adapted to the constained case, either 


by considering the N — 1 first equations in (2.4), or by solving a constrained optimisation 


problem. In particular, the stationary points have to satisfy 

W 0 (®) = A1 


(3.2) 


for a Lagrange multiplier A E M (this is indeed consistent with (2.6)). In addition, the 
constraint x € S has to be satisfied. 


In components, the condition (3.2) becomes 


- Xi = A 


i = 1, • • • ,N . 


(3.3) 


Let A c = ^ 7 . The equation £ 3 — £ = A has three real solutions if | A | < A c , two real solutions 
if | A| = A c and one real solution otherwise. The last case is incompatible with the constraint 
while the second case can only occur if N is a multiple of 3, because then the two 
solutions of £ 3 — £ = A have a (—2 : 1 ) ratio. 

We henceforth assume that |A| < A c , and denote by ao,a\,a 2 the distinct roots of 
£ 3 — £ — A. Then each 37 solving (3.3) has to be equal to one of the ctj. We let o ? - be the 


6 







Figure 3. Transition graph of the constrained system for N = 4 and 7 = 0. Black and 
white circles represent respectively coordinates equal to 1 and to —1. A few 1-saddles 
associated with edges of the graph are shown, with blue circles indicating coordinates equal 
to 0. 

number of occurrences of ctj , and reorder the a 3 in such a way that ao ^ a\ ^ ai- We 
denote such a stationary point by the triple ( 00 , 01 , 02 )- Observe that we necessarily have 
ao + 0,1 + a 2 = N. 

Proposition 3.1 (Local minima and 1-saddles for 7 = 0). Assume that N is not a multiple 
of 3, and let x* be a critical point with triple ( 00 , 01 , 02 ). Then 

• if 2 ai >00 + 02 , then x* is a stationary point of index ao; 

• if 2a\ <00 + 02 , then x* is a stationary point of index 02 — 1 . 

We give the proof in Appendix |A.1| It is based on the construction of an orthogonal 
basis around each stationary point, in which the Hessian matrix is block-diagonal with 
blocks of size 3 at most, so that the signs of its eigenvalues can be determined. 

Remark 3.2. The case 2ai = 00 + 02 can only occur if N is a multiple of 3, because 

ao + 02 = N — a\ would imply ai = N/3. In case IV is a multiple of 3, there exist 

one-parameter families of degenerate stationary points HZ3- For simplicity we exclude this 
situation in all that follows. 0 

Proposition |3.1| yields the following classification of local minima and saddles of in¬ 
dex 1 : 

1. Local minima x * £ So necessarily have triple (0, a, N — a) with N/3 < a ^ N/2. 

2. Saddles of index 1 either have triple (l,o, IV — a — 1) with N/3 < a (IV — l)/2, or 

they have triple (IV — 2 — a, a, 2) with N/2 — 1 ^ a ^ 2 and a < N/3. The latter case 

can only occur if N = 4, and corresponds to the triple (1,1, 2). 

Example 3.3 (The case N = 4). If N = 4, then So contains 6 points, consisting of 

all possible permutations of (1,1,—1,-1). In addition, there are 12 saddles of index 1, 
consisting of all possible permutations of (1, —1,0,0). Each of these saddles connects the 
two local minima obtained by replacing one 0 by 1 and the other one by — 1 , and vice 
versa Isa na Section 2.4]. The associated transition graph is an octahedron (Figure [3j. ♦ 
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Figure 4. Example of transition rules for N = 8. The coordinates for family Bq are • = 1 
and O = —1. Those for Cj are O = —1/\/7, O = 3/\/7 and O = —2 /a/ 7. Those for Bi are 
O = 5/-\/l9 and O = —3/-/T9. 


We will henceforth limit the discussion to the case where N = 2 M is even, N Js 8 and 
N is not a multiple of 3. Then the 1-saddles necessarily correspond to triples of the form 
(1, a, N — a — 1). In order to ease notations, we write fc max = 6j and 

• Bk for the set of all local minima with triple (0, M — k, M + k), where k G [0, A; max ]]; 

• Ck for the set of all 1-saddles with triple (1, M — k,M + k — 1), where k G [1, A; max ]]. 
Simple combinatorics shows that the cardinalities of these families are 


I-Bn I = 


2 M 
M 


\B k \ = 2 


2 M 
M + k 


\C k \ = 


2(2 M)\ 


1 \(M - k)\(M + k- 1)! 


(3.4) 


where k G [1, A; max J. The factors 2 are due to the fact that except for Bo, there are always 
two choices for the signs of coordinates. 

One can obtain explicit expressions for the coordinates of all these stationary points, 
see (A.3) in Appendix | A. 1| Here it will suffice to know that local minima in Bq simply have 
M coordinates equal to +1 and M coordinates equal to —1. These stationary points are 
expected, and admit a simple interpretation in terms of a particle system: we just associate 
each coordinate equal to +1 with the presence of a particle, and each coordinates equal to 
— 1 with the absence of a particle, that is, a hole. 

The other families of local minima B1 ,... , Bfc max have more complicated coordinates, 
which do not allow for an interpretation as a particle system. In fact their presence comes 
a bit as a surprise, so that we will call them spurious configurations. We will however show 
below that they have a higher energy than the configurations in Bo, and therefore they will 
not play an important role when the system is observed on a sufficiently long timescale. 


Example 3.4 (The case N = 8). If N = 8, there are two families of local minima Bq and 
B1, and one family of 1-saddles C\ (Figure [4]) . 

• The family of local minima Bo corresponds to the triple (0,4,4), and contains all points 
that have 4 coordinates equal to +1 and 4 coordinates equal to —1. They can thus be 
interpreted as configurations with 4 particles and 4 holes. 

• The family of local minima B\ corresponds to the triple (0, 3, 5). It contains all points 
with 3 coordinates equal to 5/x/l9 and 5 coordinates equal to —3/\/l9, as well as all 
configurations with opposite signs. 

• The family of 1-saddles C\ corresponds to the triple (1, 3,4). It contains all points with 
1 coordinate equal to —1/\/7, 3 coordinates equal to 3/\/7 and 4 coordinates equal to 
-2/V7, as well as all configurations with opposite signs. ♦ 


Now that we have determined all stationary points in Bo and <Si, we have to find the 
structure of the transition graph Q = (So,£). In other words, we have to determine which 
local minima are connected by a given 1-saddle. This question is answered in the following 
result. 
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Theorem 3.5 (Transition graph for 7 = 0). Each 1-saddle in C-k connects exactly one local 
minimum in Bk-i with one local minimum in Bk- More precisely, if the coordinates of 
the saddle have values Oq, a \, a' 2 , and those of the local minima are respectively 07 , 0:2 and 
then the connection rules are given by 


a\ v- 


-> «2 

1 coordinate , 


a\ v- 

— > a'i V- 

->• a" 

M — k coordinates , 

(3.5) 

a 2 <- 

— >■ a 2 

-> a 2 

M + k — 1 coordinates . 



We give the proof in Appendix |A.1| It is based on the construction of two continuous 
paths connecting a given point in Ck to one point in Bk -1 and one point in Bk , such that the 
potential decreases along the path when moving away from the saddle. Figure [4] illustrates 
the connection rule in the case N = 8. See also U Fig. 5], 


Using the relations (3.4), one easily checks that the number of saddles in Ck is indeed 


equal to the number of allowed connections between elements in Bk -1 and Ck as well as Bk 
and Ck- 


3.3 The case of weak positive coupling 

It follows from basic perturbation arguments that the transition graph Q will persist for 
small positive coupling intensity 7. Indeed, if we assume that N is not a multiple of 3, then 
all stationary points for 7 = 0 are nondegenerate, so that the implicit function theorem 
shows that they still exist for small positive coupling, and move at most by a distance of 
order 7. In addition, perturbation results for the eigenvalues of matrices such as the Bauer 
Fike theorem (see for instance [21]) show that the signature of nondegenerate stationary 
points does not change for small 7. Finally, the proof of Theorem |3.5| essentially relies on 


the relation (A. 10), whose coefficients depend continuously on 7. 


The drawback of this argument is that while it shows that for any N < 00, there 
exists a critical coupling 7* (IV) > 0 such that the transition graph does not change for 
0 ^ 7 < 7* (IV), it does not yield a good control on the critical coupling as N —> 00. To 
obtain a lower bound on 7*(IV) which is uniform in N (at least for k fixed), we adapt 
from J6] an argument based on symbolic dynamics to obtain the following result. 

Theorem 3.6 (Persistence of the transition graph for small positive 7). There exists a 
constant c > 0, independent of N, such that the stationary points of the families Bk and Ck 
persist for 

, (3.6) 

without changing their index. In the particular case of stationary points of the family Bq, 
we have the sharper result that they persist at least as long as 7 < 7 — V5 — 0.097. 


The proof is given in Appendix |A.2[ It also provides a criterion allowing to sharpen the 
bound (3.6) for families other than Bq, cf. (A.38), which however is not essential in what 
follows. 

The important aspect of this result is that all families of stationary points Bk or Ck with 
bounded away from g are ensured to exist up to a positive critical coupling independent 
of N. Only stationary points with k = ^ — o(N ) might disappear at a critical 7 which 
vanishes in the large -N limit. 
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4 Metastable hierarchy 


Now that the structure of the transition graph Q is understood, we have access to information 
on timescales of the metastable process. A convenient way of doing this relies on the concept 
of metastable hierarchy , which is an ordering of the local minima from deepest to shallowest. 


We summarise this construction in Section|4.1[ before applying it to our case in Section 4.2 


A more refined hierarchy can be obtained for small positive coupling 7 among the local 


minima of the family Bo, which have a particle interpretation; we do this in Section 4.3 


4.1 Metastable hierarchy and Eyring—Kramers law 


We consider in this section a general reversible diffusion process in 
with potential V of class C 2 . 


of the form ( 2 . 2 ), 


Definition 4.1 (Communication height). Let x* be a local minimum ofV and let A C M^. 
The communication height from x * to A is the nonnegative number 


H(x*,A)= inf sup V(jft)) — V(x*) , (4.1) 


where the infimum runs over all continuous paths 7 : [ 0 , 1 ] 


t>N 


such that 7 ( 0 ) = x* and 


7 ( 1 ) E A. Any path 7 realising (4.1) is called a minimal path from x * to A. 


The communication height measures how high one cannot avoid climbing in the potential 
landscape to go from x* to A. Assuming A does not intersect the basin of attraction of 
x* and all stationary points of V are nondegenerate, it is not difficult to show that the 


supremum in (4.1) is reached at a 1-saddle z* of V (see for instance [ 8 ], Section 2]). In that 
case, one has H(x*,A) = V{z*) — V(x*). 

A notion of metastable order of local minima was introduced in HU- In our case, due to 
the fact that many minima have the same or almost the same potential value, we introduce 
the following generalisation of this concept to partitions of the set of local minima. Typically, 
we will apply this definition to cases where the points in each element of the partition have 
approximately or exactly the same potential height. 

Definition 4.2 (Metastable hierarchy of a partition). A partition Sq = Pi U P 2 hi... U P m 
of the set So of local minima of V is said to form a metastable hierarchy if there exists a 
constant 6 > 0 such that for all k E [2, m\, one has 


Hi x 


k -1 

U p < 

2—1 


min H 

y*£P e 


*,U Pi\Pt 


2=1 


for all x * E Pk and all £ E [1, k — 1]. In this case, we write 


-6 


(4.2) 


P\ P P2 A ■■■ P P-n 


(4.3) 


In words, it is easier, starting in any point in Py, to reach a lower-lying set Pg in 
the hierarchy than it is, starting in such an Pg, to reach any other set among P\,... Pk■ 
A graphical way of constructing the hierarchy relies on the so-called disconnectivity tree 
H2] ; it is illustrated in Figure [ 5 ] in a simple case where all Pk = {x£} are singletons. 
The leaves of the tree have coordinates (x k ,V(x k )); each leaf is connected to the lowest 
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Figure 5. Example of a 4-well potential, with its disconnectivity tree. The metastable order 
is given by x\ -< x\ -< x\ -< x\. The communication heights H k = H(x k ,{x *,... ,x k _ 1 }) = 
V(z k ) — V(x k ) provide the Arrhenius exponents for mean transition times and small eigen¬ 
values of the generator C. Prefactors in the Eyring-Kramers law (4.4) are given in terms of 
second derivatives of the potential at the local minima x k and 1-saddles z k . 


saddle reachable from it, and the procedure is repeated after discarding the shallower local 
minimum whenever two branches join. 

In the particular case where all P k are singletons, the following result by Bovier, Gayrard 
and Klein connects the metastable hierarchy with certain first-hitting times and with small 
eigenvalues of the infinitesimal generator C = eA — VK(x) ■ V of the diffusion. 

Theorem 4.3 (Eyring-Kramers law for nondegenerate potentials |11]). Assume the local 
minima of V admit a metastable order x* -< ■ ■ ■ -< x^. For each fee [1, m], denote by r k 
the first-hitting time of the e-neighbourhood of {x^,..., x* k }, and let X k by the kth smallest 
eigenvalue of —C. Assume further that for each k, there is a unique 1-saddle z k such that 
any minimal path from x£ to { x\ ,..., xt 1 } reaches communication height only at zt- Then 
for each k e [[2, ttt,]] ^ one has 


® Xk [T k -l] = 


2ir 


I A- (z* k ) | 


/detvwll' e '' /W> —'^‘l 1 + 0(s ,/2 Hog E |^)] , (4.4) 


where V 2 K(x) denotes the Hessian matrix of V at x, and A ~(z k ) is the unique negative 
eigenvalue of\7 2 V(z k ). Furthermore, X\ = 0 and there exists a constant 9\ > 0 such that 


Afc — 


e 


[1 + 0(e- ei / £ )] 


(4.5) 


holds for all k e [2, m\. 

This result tells us in particular that if the system starts at a stationary point at the 
end of the metastable hierarchy, it will spend longer and longer amounts of time going 
down the hierarchy (possibly visiting other local minima in between), before reaching the 
ground state x*. In particular, the spectral gap A 2 — Ai = A 2 of the system, which gives the 
exponential rate of convergence to equilibrium, depends to leading order only on the second 
local minimum in the hierarchy x|, and on the saddle z\ connecting it to the ground state. 
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Figure 6. Value of the potential V 0 along a path B 0 —>■ C\ — > B 1 — > ... in the case 
N = 20 (not to scale). The associated disconnectivity tree shows that the Bk are indeed in 
metastable order. Thus the long-time dynamics will concentrate on the set Bq of particle- 
hole configurations. 


4.2 Hierarchy on the families Bk 


Unfortunately, Theorem T3 does not apply to our situation, because one cannot find a 
hierarchy for singletons. This is due to the fact that the potential Uy has many symmetries, 
and therefore many stationary points have the same potential height, preventing us from 
fulfilling (4.2) with a positive 6. In particular, in the uncoupled case 7 = 0, the system 


is invariant under the group G = & y x Z 2 , where is the symmetric group describing 
permutations of the N coordinates, and the factor Z 2 = Z/2Z accounts for the x i —> —x 
symmetry. The families Bk and Ck each form a group orbit under G , that is, they are 
equivalence classes of the form {gx : g G G}. 


However, we will be able to draw on results of [5], which generalise Theorem 4.3 to 
Markovian jump processes invariant under a group of symmetries, and the extension of 
these results to diffusion processes mm- In particular, [5j Thm 3.2] shows that if the 
system starts with an initial distribution which is uniform on some Bk, then a very similar 


result to Theorem 4.3 holds true. The only difference is that the prefactor in the Eyring- 


Kramers law (4.4) has to be multiplied by a factor which can be explicitly computed in 


terms of stabilisers of the group orbits. 

The following result provides a metastable order on the Bk, which is exactly what is 
required to apply the theory from 13 EZl EE! in the uncoupled case. 

Theorem 4.4 (Metastable hierarchy on the Bk). If 7 = 0, then the families Bk satisfy a 
metastable order given by 

Bo -< Bi -< ■ ■ ■ -< Hfc max . (4.6) 

Furthermore, any minimal path from Bk to Bk -1 reaches communication height only on 


saddles in Ck- The hierarchy (4.6) still applies for sufficiently small positive 7 , the only 
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difference being that all points inside a given Bk do not necessarily have the same potential 
value. 


We give the proof in Appendix B.l The situation is illustrated in Figure [ 6 } As k 
increases from 0 to & max , the potential height of the Bk increases, while the barrier height 
between Ck and Bk decreases. See Appendix |B.1| for explicit expressions for these potential 
values. Applying [5j Thm 3.2], we obtain in particular the following result. 


Corollary 4.5. Fork £ [l,A; max ]], iy ._1 be the first-hitting time of the £-neighbourhood of 
Bo U • • • U Bk~ i. If the initial distribution p of the system is concentrated on Bk U • • • U -Bfc max 
and invariant under G, then for 7 = 0 one has 




2vr 


\X-(4)\(M + k) 


\j det wT°K) [1 + O(e 1/2 |log T /2 )] . (4.7) 


where x\ is any local minimum in Bk, z k is any saddle in Ck, and M = N/2. 

PROOF: Theorem 3.2 in [5] shows that in the case of a symmetric initial distribution, the 

has to be multiplied by the factor | G x * nG x *_J/|G x *|, 
where G x = {g £ G : g(x) = x } is the stabiliser of x. If k 1, then \G x * k \ is the number 
of permutations that leave invariant any element in Bk, and is equal to (M — k)\(M + k)\. 
Similarly, | G x * n G x * | is the number of permutations leaving invariant any two elements in 
Bk and Bk -1 connected in the transition graph Q, which is equal to (M — k)\(M+k — 1)!. □ 


usual Eyring-Kramers formula (4.4 


Note the extra factor (M + A ;) -1 in (4.7). In fact, M + k is also the number of saddles in 
Ck that are connected with any given element of Bk (cf. {5J Eq. (2.25)]). The interpretation 
of this factor is that since the system has M + k different ways to make a transition from 
a given x* k £ Bk to Bk-i, the transition time is divided by this factor. 

The above result will still apply for small positive coupling, but with a more complicated 
expression for the prefactor. This is because the system is no longer invariant under XZ 2 , 
but under the smaller group 0 jv x Z 2 , where Dn is the dihedral group of symmetries of 
a regular Wgon. The important aspect for us is that we still have a control of the time 
needed to reach the family of stationary points Bo, which lie at the bottom of the hierarchy 
and have an interpretation in terms of particle-hole configurations. The dynamics among 
configurations in Bo is much slower than the relaxation towards Bo, because it involves 
crossing the potential barrier from Bo to B\ via C\ . We will analyse it in more detail in 
the next section. 


4.3 Hierarchy on Bq and particle interpretation 

We assume in this section that 0 < 7 < 7 C , where 7 C is the critical coupling below which 
all stationary points in Bo, B\ and C± exist without bifurcating. The central observation 
in order to classify points in Bo is that if x*(y) is any critical point of V^, then 

N 

E 7 (.t*( 7 )) = Fo(s*(0)) + J^K +1 (°) - x* (0 )) 2 + 0 ( 7 2 ) ■ (4.8) 

i— 1 


This is because Vo(®*( 7 )) = Vo(a;*(0)) + G( 7 2 ), as the first-order term in 7 vanishes since 
VVo(x*(0)) = A1 is orthogonal to x*(j)—x*(0), which belongs to the hyperplane S. The first 
term on the right-hand side of (4.8) is constant on each Bk and each Ck- The second term 
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Figure 7. Example of an allowed transition, from a configuration in B 0 with two interfaces 
to a configuration in Bq with 4 interfaces. The net effect is that a particle has hopped by 
two sites. 

is determined by the number of nearest-neighbour coordinates of x*(0) that are different, 
which we are going to call interfaces of the configuration. 

In particular, if x*(0) £ Bq, we know that all its components have values ±1. We define 
its number of interfaces as 


N 



(4.9) 


i=1 


so that we have 



Interpreting each 1 as a particle and each —1 as a hole, it is natural to introduce the 
following terminology: 

• a site i with 2 interfaces will be called an isolated particle or hole; 

• a sequence of at least 2 contiguous particles or holes will be called a cluster-, 

• a site with 1 interface lies at the boundary of a cluster ; 

• a site without interface belongs to the bidk of a cluster. 

Lemma 4.6. Let x * be a critical point in Bq and write M = ^ ^ 4. Then the following 
properties hold. 

1. The total number of interfaces is even. 

2. If I\/_\{x*) = 2, then x* consists in a cluster of M particles and a cluster of M holes. 

3. Ifly-xix*) > M, then x* has at least one isolated site. 

4. Among the x* £ Bq with J 1 /_ 1 (x*) £ [4, M\. there exist both configurations with isolated 
sites and configurations without isolated sites. 

PROOF: Denote by N c the number of clusters, by N, the number of isolated sites, and by 
p = hz-iix*) the number of interfaces. Then we have p = N c + N,, which is necessarily 
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even. Since clusters have at least two sites, N ^ 2N C + N{, implying N c ^ N — p and thus 
N\ ^ 2p — N. Thus if p > M, then N[ > 0. If 4 ^ p ^ M. then a possible configuration 
consists in p — 2 clusters of size 2, leaving at least 4 sites that can be split into 2 more 
clusters. Another possibility is to have p— 2 isolated sites, leaving at least N — 2 sites that 
can again be split into 2 clusters. If p = 2, we necessarily have 2 clusters of equal size. □ 


This result motivates the following notation for configurations in Bq: 

• A 2 denotes the set of all configurations with interface number I\/-i(x*) = 2 ; 

• for even p E [4, M], A p denotes the set of all configurations x* E Bq with p inter¬ 
faces having at least one isolated site, and A' p denotes the set of configurations with p 
interfaces having no isolated site; 

• for even p E \M + 1, NJ, A p denotes the set of configurations with p interfaces (which 
all have at least one isolated site). 


We now need to determine the communication heights between configurations in these 
different sets for small positive 7 . For this, we have to take into account the fact that any 
transition between two configurations in Bq involves crossing two 1-saddles in C \, separated 
by an element of B\ (Figure [7J). The communication height will thus be determined by the 
highest of the two saddles. Examining the different possible cases yields the following result, 
which is proved in Appendix B.2| 


Proposition 4.7 (Transitions between configurations in Bq). Let x\fiy), x\fy) E Bq be 
two particle/hole configurations, and denote by p = Ii/_i(xi( 0 )) the number of interfaces 
of x\{ 0). Then a transition between these configurations is possible if and only if x%( 0) is 
obtained by interchanging a particle and a hole in x*(0). The interface number of xfiO) 
satisfies 

li/_i(x£(0)) E {p-4,p-2,p,p + 2,p + 4:} . (4.12) 

The communication height from x\{pf) to x^fiy) admits the expansion 


H { xU 7),*2(7)) = # (0) +7tf (1) (*i(0),s£(0)) +0(7 2 ) , (4-13) 


where 

H { ° ] = ~ 4(m( ( -3M 1 !3) (414> 

depends only on M = while H^^x*^), x^O)) also depends on p and on the number of 
interfaces of the two exchanged sites as detailed in Table [7} 


Table [l] shows that all allowed transitions between particle/hole configurations have 
simple physical interpretations. In particular, only the last four types of transitions decrease 
the number of interfaces. Types V.b and V.c can be viewed as an isolated particle merging 
with another particle (isolated or at the boundary of a cluster), type V.a as a particle 
splitting from another one to fill a hole between two particles, and type VI as an isolated 
particle jumping into a hole between two particles. Types I and II are just the reversed 
versions of types VI and V, while all transitions of type III and IV are their own reverse. 

Figure [ 8 ] shows the allowed transitions in the case N = 8 ; only transitions that minimise 
the communication height are shown. Figure [9] shows the case N = 16. Note that in 
accordance with Lemma 4.6, only configurations with p M interfaces appear in the two 
types A p (with isolated particles and/or holes) and A' p (without isolated particles and/or 
holes). 
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Transition 

Ap 

# (1) (*i( 0),x£(0)) 

Saddle 

I 

••• 

+4 

10AL 2 - 36 M + 36-3 p 

[0,2,p + 2] 

4(M 2 - 3 M + 3) 

II.a 

••• 000^^060 ■■■ 

+2 

2(M - 3) 2 - 3 p 

4(A/ 2 - 3M + 3) 

[0,2 ,p] 

Il.b 

00 O ■■■ OOO 




II.c 

A 

... 0000 ••• 




III 

••• 

0 

-2 M 2 + 6 M - 3 p 

4(M 2 - 3 M + 3) 

[1, l,p - 1] 

IV.a 

••• ••• 

0 

-6M 2 + 12 M - 3 p 

4(A/ 2 - 3 M + 3) 

[0,2,p-2] 

IV.b 

••• o»o ••• OOO ••• 




IV.c 

A 

... 0000 ■■■ 




IV.d 

A 

• •• 0000 ••• 




V.a 

••• ••o^-~'#o# ••• 

-2 



V.b 

• ■ ■ 000 • • • 'oo» • • • 




V.c 

A 

... o«0# ••• 




VI 

••• O 0 O^~^ 0 O 0 ■■■ 

-4 




Table 1. List of allowed transitions between elements of Bq , viewed as a particle moving 
into a hole. The different columns show, respectively, the type of transition, the change A p 
of the number of interfaces, the first-order correction to the communication height, and the 
numbers of interfaces of types a^/a^, ot^/a^ and c^/c^ of the highest saddle encountered 
during the transition (cf. Appendix |B.2[). 


The first-order correction H 4) to communication heights depends not only on the num¬ 
ber M of particles, but also on the number p of interfaces. This is a nonlocal effect of the 
mass-conservation constraint. However, in the limit M —> oo, the four possible corrections 
converge respectively to §, |, —\ and — i.e. they no longer depend on p. 

With this information at hand, it is now possible to determine the metastable hierarchy 
among the families A p and A' p . The result, which is proved in Appendix 
follows. 


B.2 


reads as 


Theorem 4.8 (Metastable hierarchy of particle/hole configurations). Let M' be the largest 
even number less or equal M = y. Then 

A 2 A A 4 + Aq -<■■■-< Am'_ 2 A A'mi + A 4 -< Ag + • • • + Ajy —2 A Ajy . (4.15) 

defines a metastable order of the families A p and A' p . 
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6 III 



Figure 8. Minimal transitions between particle/hole configurations in B 0 for N = 8. Ar¬ 
rows indicate transitions that decrease the energy, and are labelled according to Table [T| 
Each node displays only one representative of an orbit for the group action of T>n x Z 2 . The 
other elements of an orbit are obtained by applying rotations, reflections and interchanging 
particles and holes. Blue nodes represent stationary points in B\. Not shown are transitions 
within the families A p and A' p , which are of type III or IV. 


5 Analysis of the dynamics 


5.1 Interface dynamics 

The transition rules and communication heights given in Proposition [4V] and the metastable 
hierarchy obtained in Theorem |4.8| yield complementary information on the dynamics be¬ 
tween particle/hole configurations in Bq. Recall that the process behaves essentially as 
a Markovian jump process with transition rates of order , while the hierar¬ 


chy (4.15) classifies the states according to the time the process spends in them in metastable 


equilibrium. 

At the bottom of the metastable hierarchy, we find the set A 2 of configurations having 
one cluster of M particles: this constitutes the ground state of the system, which can be 
interpreted as a solid or condensed phase. At the top of the hierarchy on Bq, we find the set 
An of states with N interfaces. These consist in M isolated particles, and can be interpreted 
as a gaseous phase. 

The transition graph implied by Proposition 4.7 (and illustrated in Figures [8] and [9]) 
shows that when starting in the configuration An, the most likely transitions gradually 
decrease the number p of interfaces, in steps of 2 or 4. Thus the system tends to gradually 
build clusters of increasing size. As the communication heights given in Table [l] increase as 
p decreases, this condensation process becomes slower as the size of clusters increases. This 
is different from the usual Kawasaki dynamics, in which the transition rates depend only on 
the change Ap of the number of interfaces. Note in particular that for given p, transitions 
of type IV, V and VI all occur at the same rate. 
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Figure 9. Minimal transitions between particle/hole configurations in Bo in the case N = 
16. Arrows indicate transitions that decrease the energy. 


When the number p of interfaces reaches M (meaning that there are on average 2 par¬ 
ticles per cluster), new configurations A' p become possible. These consist of p clusters sepa¬ 
rated by at least 2 sites, and appear as dead ends on the transition graph. The metastable 


order (4.15) shows that these configurations are actually more stable than those of type A p , 
p ^ 4, which have isolated particles or holes, and act as gateways to configurations with 
fewer interfaces. In particular, configurations in A' a are those with the longest metastable 
lifetime. The system can spend considerable time trapped in configurations with p ^ 2 
clusters of particles, separated by p clusters of holes (as seen in Figure [l]). 


5.2 Spectral gap 

Another interesting information on the process that can be obtained from its metastable 
hierarchy is its spectral gap. We already know that the generator C admits the eigenvalue 0, 


which is associated with the invariant distribution (2.9). This eigenvalue is simple because 
the process is irreducible and positive recurrent. The spectral gap is thus given by the 
smallest nonzero eigenvalue A 2 of — £, which governs the rate of relaxation to equilibrium. 

At first glance, one might think that the spectral gap has order e ^ 7(2/*))/^ where 

z* is a 1-saddle in C\ and y* is a local minimum in B\. Indeed, this is the inverse of the 
longest transition time obtained in Corollary 4.5. However, the corollary only applies to 


symmetric initial distributions, and transitions from B\ to Bq via C\ are not the slowest 
processes of the system. In fact, this role is played by transitions between configurations 
in A 2 , which occur via saddles in C 1 , leading to a spectral gap of order e~^ ( z *l~ v r ( x *))/- ; 
where x* is a local minimum in A 2 rather than B\. Applying the theory for symmetric 
processes in [5], we obtain the following result. Its proof is given in Appendix [C| 


Theorem 5.1 (Spectral gap). If £ is small enough, then the smallest nonzero eigenvalue 
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of —C is given by 


A 2 = 4 sin 2 — 


\-U 


' det V 2 i4y(x*) 

|det V 2 U 7 (.s*)| 


e -[^(**)-^(**)]/ C [ 1 + 0 ( £ l / 2 | log £ | 3 / 2 jj ( ( 57 ) 


where x * is any configuration in A 2 , and z* is any saddle in C\ whose limit as 7 
exactly 3 interfaces. In particular, we have 

T .,*v M(M - 1) M 2 — 6M + 6 . 2 . 

7 ^ ) ) “ 4(M 2 -3M + 3) + 7 2(M 2 - 3M + 3) + ^ 

= | + ^7 + 0 (iV' 1 ) + 0(7 2 ). 


Furthermore, 


/ detV 2 U 7 (x*) /- 
|A_( ~ |detV 2 F 7 ( 2 *)| V2 


M 2 - 3M + 3 


(M- l)^/M{M- 3) 
= >/2 + ^(iV- 1 ) + 0 ( 7 ) . 


M—2 


+ 0 ( 7 ) 


0 has 


(5.2) 


(5.3) 


The fact that the spectral gap (|5.lh decays like IV 2 for large N is highly nontrivial. It 


is related to the fact that the symmetry group T>n x Z 2 admits irreducible representations 
of dimension 2 , and its computation requires the full power of the theory developed in [5], 
Physically, this result means that some transitions between states in A 2 require a time 
of order N 2 e 1//4e , i.e., increasing as the square of the system size when the noise intensity e 
is constant. In other words, the motion of interfaces slows down like N~ 2 when the system 
becomes large. 


6 Conclusion 


Let us briefly summarise the main results obtained in this work. 

• Using the concept of metastable hierarchy, the long-term dynamics of the system can 
be reduced to an effective process jumping between particle/hole configurations. These 
configurations exist as long a the coupling intensity 7 is smaller than a critical value, 
bounded below by a constant independent of the system size. 

• The effective dynamics tends to reduce the number of interfaces, and slows down as 
this number decreases. As soon as the average size of clusters reaches 2 , the system 
can get trapped in configurations without isolated sites, which are more stable than any 
configuration with isolated sites. 

• The spectral gap is of order IV -2 e' -1 / 4 ^, which decreases as the square of the inverse of 
the system size. This means that transitions between the N configurations forming the 
ground state A 2 slow down as N increases. 


We emphasise that all results obtained here apply for arbitrarily large but finite system 
size N. In fact, some quantities like the number 8 defining the metastable hierarchy go to 
zero in the limit N —» 00 , so that the orders ( |4.6[ ) and ( |4.15[ ) only make sense for finite N. We 
do not claim either that the error terms of order e 1 / 2 |loge | 3 / 2 in (5.1) and (4.7) are uniform 
in N, though results obtained in a similar situation in [3] indicate that they probably are. 

A different situation of interest, not considered here, arises when the coupling intensity 
7 grows like N 2 . Then one expects that the system converges to a mass-conserving Allen- 
Cahn SPDE on a bounded interval, which has considerably fewer metastable states. Indeed, 
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an analogous scenario was obtained in [?j, where the unconstrained system with 7 ~ N 2 
was shown to have only 2 local minima, and at most 2N saddles of index 1. If, by contrast, 
one has 1 <C 7 <S N 2 , a scaling argument shows that the system should converge to 
an Allen-Cahn SPDE on a growing domain, which admits more metastable states; see in 
particular EDES] for results in the unconstrained case, and pB] for a recent convergence 
result in dimension 2 . 

The behaviour of the constrained system for lattices of dimension larger than 1 remains 
so far an open problem. The phenomenology is expected to be different, because the energy 
of clusters then depends not only on the size of their interfaces, but also on the size of their 
bulk. This can result in scenarios where the interface dynamics accelerates once a critical 
droplet size has been reached, as is well known for lattice systems with standard Kawasaki 
dynamics [14] . 


A Proofs: Potential landscape 


A.l The uncoupled case 


Proof of Proposition |3l7| Consider a critical point x * of the constrained system with triple 
(ao, ai, CI 2 ). Recall that this means that x* has a,j coordinates equal to 07 , j = 0,1, 2, where 
the Qj are distinct roots of £ 3 — £ — A for some A £ (—A C ,A C ). By Vieta’s formula, these 
roots satisfy 

ao T a 1 T a 2 = 0 . (A.l) 


We always have ao + ai + 0,2 = N, and by convention ao ^ a\ ^ 02 - Note that we may 
assume ao / 02 , since otherwise all aj would be equal, and thus N would be a multiple of 
3, which is excluded by assumption. 

Combining (A.l) with the constraint ^ x* = ao<so + aiai + 0207 = 0 yields the relation 


(ai — ao)aq + (02 — 00)^2 — 0 • 


(A.2) 


Solving for <37 and using the fact that all a 3 — 07 are equal, a short computation shows that 

ao = ±(ai — , 

a \ = ±(02 — ao )-?? 1 / 2 , ( A . 3 ) 

a 2 = ±(«o — ai )-?? 1 / 2 , 


where 


R = 


®2 T ai — 2ao 
(ai - a 0 ) 3 + (a 2 - a 0 ) 3 


1 

a o + a i + a 2 _ a o a i — ao°2 — a\Q2 


(A.4) 


We now turn to determining the signature of the Hessian at these critical points of the po¬ 
tential Ky restricted to the hyperplane S. This signature does not depend on the parametri- 
sation of S , so that it is equal to the signature of the Hessian of 


R 7 0i, • • •, xn-\) = Vy ( Xl ,..., xn-i, —xi — ■ ■ ■ — xn- 1 ) 
Computing the Hessian of Ky at x* shows that it has the form 

0 0 \ /! 


'(3a 2 - l)ll ao 
H= | 0 

0 


0 

(3a 2 — l)l ai 0 \+(3a 2 -l) 

0 (3a 2 — l)lla 2 -i / 


V 1 


(A.5) 


(A. 6 ) 
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where l a denotes the identity matrix of size a. We now distinguish between the following 
cases. 


1. 

2 . 


( 00 , 01 , 02 ) = (0,0,/V). Then x* = 0, and one easily sees that —H is positive definite, 
so that x* is a saddle of index N — 1. 

ao = 0 and ai ^ 1. Using the expressions (A.3), we obtain that 3af — 1 > 0 and 
3«2 — 1 has the same sign as 2oi — 02 - Let ei,... ,e^-i denote the canonical basis 
vectors. Then {ei — are eigenvectors of H with eigenvalue 3af — 1, and 

{e ai +i — [a 1 + 2 ,JV— 1 ] are eigenvectors of H with eigenvalue 3 a\ — 1. 

To hnd the remaining two eigenvalues, let u = Yli =1 e * an d v = X^ai+i e «- These two 
vectors span an invariant subspace of H, in which the action of H takes the form 


Ai 


f 3 af — 1 + ai(3a2 — 1) (02 — l)(3o| — 1)\ 
y ai( 3 a| — 1 ) 02(302 — 1 ) ) 


(A.7) 


Computing the determinant and the trace of Ai, one sees that if 2ai > 02 , then the 
two eigenvalues of Ai are strictly positive, so that x* is a stationary point of index 0. 
If 2ai < 02 , then Ai has one strictly positive and one strictly negative eigenvalue, and 
x* has index 02 — 1. 

3. o 0 ^ 1. In that case one finds that 3of — 1 > 0, while 3 oq — 1 has the same sign as 
(2o2 —ai—oo)(ao — 201 + 02 ) and 3a 2 — 1 has the same sign as (2ao—ai — O 2 )(ao — 201 + 02 ). 
Here it is better to invert the roles of a\ and 02 in the expression for H . Similarly to the 
previous case, one finds ao — 1 eigenvectors with eigenvalue 3— 1, 02 — 1 eigenvectors 
with eigenvalue 3a 2 —l and oi— 2 eigenvectors with eigenvalue 3af — 1 (these eigenvectors 
are of the form ei — e*, e ao +i — e* and e ao+a2 + 1 — e * f° r appropriate ranges of i). 

To hnd the other eigenvalues, let u = YaI 1 e;, v = 1 e * and w = Z)iIao+a 2 +i e *- 

These span an //-invariant subspace, in which the action of H takes the form 

/3«q — 1 + oo(3af — 1) a2(3af — 1) (ai — l)(3a^ — 1)\ 

M. = I ao(3af — 1) 3«2 — 1 + 02(3af — 1) (ai — l){3a\ — 1) I . (A.8) 

\ ao(3a^ — 1) a2(3a^ — 1) ai(3af — 1) / 


In this case, one finds TV M. > 0, and det M. has the same sign as a 0 — 2ai + 02 - If 
det A4 < 0, then Ai has two strictly positive and one strictly negative eigenvalue, and 
x* has index ao- If det Ai > 0, computing the term of degree 1 of the characteristic 
polynomial of Ai one concludes that all eigenvalues of Ai are strictly positive, and that 
x* has index 02 — 1. □ 


Proof of Theorem 3.5. 

N 


Let z* 

N 


where a = y- Hl£[f + l — k max 


£ Cj, be a 1-saddle. Its triple can be written (1, a — 1, N — a) 
, yJ- We shall construct a path T, connecting z* to 


a point x* 6 B^_ 1 of triple (0, a,N — a), and such that the potential Vo is decreasing along 
r. An analogous construction holds for the connection from z * to a local minimum in B 
In fact it will turn out to be sufficient to use a linear path. Reordering the components 
if necessary, we may assume that x* = («i,..., ai, 0 : 2 ,..., 0 : 2 ) with ai repeated a times 
and a.2 repeated N — a times, and z* = (a' 0 , a\ ,..., ay, a' 2 . ■ ■ ■ ■ a^), with a\ repeated a — 1 
times and a ' 2 repeated N — a times. Note that these points indeed satisfy the connection 
rules (3.5). Let T(t) = tz * + (1 — t)x* and set h(t) = Vb(T(t)). Then a direct computation 
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shows that 


h'(t) = (a 0 - ai) [((1 - t)ai + tag) 3 - ((1 - t)a i + to' 0 )] 

+ (a — l)(a / 1 — aq) [((1 — t)a± + toy) 3 — ((1 — t)a i + iay)] 

+ (N — a)(a' 2 - a 2 ) [((1 - t)a 2 + ta' 2 ) 3 — ((1 - t)a 2 + ta 2 )] . (A.9) 

The properties of the ay and a'- yield h'( 0) = /i'(l) = 0. Since h'(t) is a polynomial of 
degree 3, it can be written as 


h'(t) = Kt(t — l)(t — i/j) 


(A.10) 


for some K,ip G M. Computing the coefficient of t 3 in (A.9) yields K > 0. Thus if we 
manage to show that -0 > 1, we can indeed conclude that h'(t) < 0 on (0,1), showing that 
h(t) is decreasing as required. The condition ^ > 1 is equivalent to having h"(l) < 0. Using 
the expressions (A.3) of the ctj, one obtains after some algebra that 

h"{ 1) = 2(u/) 2 [(a - l)(9a - 8 N) + aN 2 - a 2 N ] - AuJN(N - a )(a - 2) - u 2 aN(N - a) 
+ 3(W) 2 (iV - 0 ) [aN 3 - 3 a 2 N 2 + 3Na 2 - (a - l)(9a 2 - 9aiV + 4iV 2 )] , (A.ll) 

where co = ( N 2 — 3aN + 3a 2 )~ 1//2 and uj' = (TV 2 — 3 aN + 3(a 2 — 0 + l)) -1 / 2 stem from the 
terms i? 1 / 2 in (A.3). Using the fact that ojoj'N(N — a)(a — 2) > 0, rearranging and replacing 
oj and uj' by their values, the condition h" (1) < 0 can be seen to be true if the condition 
g(a) < 0 holds, where 

g{a) = 9(8 N - 27)o 4 + -3(56iV 2 - 156A^ - 81)a 3 + 3N(48iV 2 - 101N - 156)a 2 

- N 2 (56N 2 - 74N - 303)a + 2iV 3 (4iV 2 - 37) . (A.12) 


To check the condition, first observe that if N ^ 4 then g^\a) > 0 for all a. Next check 
that g^i^) < 0 for A ^ 4 to conclude that g^(a) < 0 for all a ^ y. Proceeding in 
a similar way with the second and first derivatives of g, one reaches the conclusion that 
g(a ) is decreasing for a ^ ^ if N ^ 4. It thus remains to show that g is negative at the 
left boundary of its domain of definition. This follows by checking the slightly stronger 
condition s(y + |) < 0 . □ 


A.2 The case of small positive coupling 

To prove Theorem |3.6j we proceed in two steps. First we ignore the constraint that station¬ 
ary points x* should belong to the hyperplane S, and prove that the equation 


VU 7 (x) = Al (A.13) 

admits exactly 3^ solutions for all ( 7 , A) in a given domain. Then we obtain conditions on 
( 7 , A) guaranteeing that these stationary points belong to S. 

Let A c = anc [ define 


D = {( 7 , A) G [0, \] x [-A c , A c ]: |A| + 76 (A) < A c (l - 7 ) 3/2 } , 


(A.14) 


where d(A) is the largest root of x 3 — x — |A|. The set D is shown in Figure 10 A simpler 
sufficient condition for being in D is obtained by observing that 


D D D’ = {( 7 , A) G [0, §] x [-Ac, Ac]: |A| < A c (l - § 7 )} , 
owing to the fact that 6(A) G [1, ^] for |A| ^ A c . 


(A.15) 
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Figure 10. The domain D in the ( 7 , A)-plane defined in (A.14) (the boundaries of D are 
not straight line segments, although they look straight). For all ( 7 , A) £ D , the equa¬ 
tion Vh 7 (i) = A1 admits 3^ stationary points. The smaller domain corresponds to the 
parameter values where stationary points of the family Bq can exist in the hyperplane S. 


Proposition A.l. If ( 7 , A) £ D, then (A.13) admits exactly 3 N 
tinuously on 7 and A. 


solutions, depending con- 


PROOF: The proof, in the spirit of |24) . is based on the construction of a horseshoe-type 
map admitting an invariant Cantor set on which the dynamics is conjugated to the full shift 
on 3 symbols. First note that we may assume 0 < 7 ^ the case 7 = 0 having already 
been dealt with. Let f\(x) = x — x 3 + A and consider the map T : M 2 —> M 2 given by 

T(x,y) = (2x-y ~-fx(x),x^ . (A.16) 

This is an invertible map, with inverse T _1 = II o T o n where II is the involution given by 
II (x,y) = (y,x). Furthermore, the relation T(x n ,x n - 1 ) = ( x n+ \,x n ) is equivalent to 


x n x n ^ (xn +1 2,x n T X n — 1 ) — A 


(A.17) 


This shows that fixed points of T N are in one-to-one correspondence with solutions of (A. 13). 
Our aim is thus to show that when ( 7 , A) £ D, the map T has exactly 3^ periodic orbits 
of (not necessarily minimal) period N. To this end, we construct some subsets of M 2 which 
behave nicely under the map T. 

We can write T(x , y) = (g(x) — y, x) where g is the function 


g(x) = 2x- -f\(x) = - [x 3 _ (t _ j)x - A] . 

'y sy L J 


(A.18) 


It has a local minimum at zq = ^/(l — 7 )/3 and a local maximum at — zq. Furthermore, it is 
strictly increasing on (—00, — zq ) and (zo,oo) and strictly decreasing on (— zq , zo ). Let a m j n 
and a max be the smallest and largest roots of x 3 — x — A. Note that max{a max , — « min } = a 
and that 

fK^bnax) — 2a max , ^(c^min) — j n . (A.19) 
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Figure 11. The sets V a and "Ho- constructed in the proof of Proposition [AT] The square is 
the set [a m in, a max ] 2 . The V a are bounded below by g(x) — a max and above by g(x) — a m i n . 
Each V a is mapped by T to the corresponding Iterating T forward and backward in 
time produces an invariant Cantor set contained in the intersections of the V a and T-L a . 


Furthermore one can check that 


( 7 , \) (E D g{~zo) > 2a max and g(z 0 ) < 2a min . (A.20) 

Denote by gZ 1 the inverse of g with range [a m i n , — zq] and introduce the “vertical” strip 

V- = {(s,y): gZ 1 (y + a min) g_ (y + ci max ), cr m m y ci max | (A.21) 

(see Figure [IT]). Then we see that T maps V_ to the “horizontal” strip T-L- = IIV-. Similarly, 
if g^ 1 denotes the inverse of g with range \—zq, zq[, then the strip 

Vo — { ( X , y ) • ^0 iy ~b ®max) X 3 q {y T Clmin)) ^min ^ 2/^5 ®max} (A.22) 


is mapped by T to Tio = II Vo- In the same way, one can construct a strip V + defined via the 
inverse gj 1 of g with range [zq, a w x ], which is mapped to T~L+ = IIV+. The property (A.20) 
ensures that the strips V CT have disjoint interiors, and the same holds for the Up. 

Consider now any finite word to = (cu_ n ,... ,to n +i) E {—,0, +} 2 ( n+1 ), and associate with 
it the set 


71+1 


I u = 


n Tk{ y -*) • 


(A.23) 


k=—n 


The above properties of the strips imply that all I u are non-empty, and have pairwise 
disjoint interior. In fact, the union of all I u converges as n —> oo to a Cantor set invariant 
under T. By a standard argument [23], for every doubly infinite sequence to E {— ,0,+} z , 
there exists an I u E [a m i n , a max ] 2 whose orbit visits V Un at time — n and 'H Wn at time n + 1 
for each n E Nq. In particular, for any of the 3^ possible ./V-periodic sequences to, we obtain 


exactly one IV-periodic orbit of T, which corresponds to one solution of (A. 13). It depends 
continuously on the parameters 7 and A, because the I u depend continuously on them. □ 

Let us point out that the above result is consistent with the previously obtained prop¬ 
erties of the system for 7 = 0 . Indeed, as 7 —> 0, the function g defined in (A. 18) becomes 
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singular, switching between —oo and +oo at the roots of x 3 — x — A, which are precisely 
the ctj introduced in Section A.l As a consequence, the invariant Cantor set collapses on 


{« 0 ) Oil, CK 2 } 2 , and the stationary points are all iV-tuples with these coordinates (there are 
indeed 3^ of them). 

In order to deal with the constraint x * E S, we will need some control on the size of the 
sets V a n T~La' • The following lemma provides upper bounds on the widths of the V CT (and 
thus also on the heights of the 'H r ?') which will be sufficient for this purpose. 


Lemma A. 2 . Assume that (7, A) E D, and denote by a m i n < a c < a max the three roots of 
x 3 — x — A. Then 

V'— d [oimiro C^min + Vl] x [« min 5 ^max] ? 

Vo C [a c - y/j, ot c + y/d ] x [a min? ^max] ? 

V+ C [a max Vt, ^max] X [c^mini ®max] • (A.24) 

PROOF: Denote by xi the x-coordinate of the top-right corner of V- (see Figure [IT]). Then 
we have the relations 




xf - (1 - y)xi - A = 7 a max , 

a (1 7)®min A — 7®min • 


(A.25) 


Taking the difference of the two lines, writing x\ = a m j n + Ai and recalling the definition 
A) = vA 1 - 7)/3 yields 


Al/ll(Ai) — 7(®max Clmin) , hi(A) — 3 (a m in A)) 4“ 3Q m inA + A" 


(A.26) 


One easily checks that the map A 1-7 /ii(A)/A is decreasing. Since x\ ^ —zq and thus 
Ai ^ —zq — a m ; n , it follows that 


^l(Ai) hi( z 0 «min) n 

A ^ — &ZQ C^min • 

Zq C^min 

As a consequence, (2zo — a m i n )Aj ^ 7(a m ax — a m i n ), so that we conclude that 

7(®max C^min) A ^ 


(A.27) 


V_ c 


C^min? C^min + 


x \a 


: mm? ^max 


(A.28) 


2^0 C^min 

Now we claim that a max ^ 2^o holds for all ( 7 , A) E D. Indeed, if g is the function defined 
in (A.18), then we have by (A. 14) 


g{ 2 z 0 ) = - [A c (l - 7 ) 3 / 2 - A] > 2d(A) 


(A.29) 


for all ( 7 , A) E D. Hence by (A. 19) we get < 7 ( 2 ^ 0 ) ^ 2a m ax = 3(a m ax)i showing as claimed 
that Q ma x ^ 2^o since g is increasing on [zq, 00 ). Using this bound in (A.28) yields the first 
relation in (A.24). 

In a similar way, if X 2 denotes the x-coordinate of the top-left corner of Vo, one obtains 
that A 2 = a c — X 2 satisfies 


A 2 MA 2 ) = 7(«max - Oi c ) , h 2 ( A) = 3 (zq - a 2 ) + 3a c A - A 2 


(A.30) 
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One obtains again that A H > / 12 (A)/A is decreasing, and its smallest value, reached at 
A 2 = a c + zo, is equal to 2zo — a c . The other relevant coordinates can be computed in the 
same way, yielding 


V 0 C 


V+ c 


OLr 


'T^max U c )\ ^ ( / y(o^c ^min) 




, OLr 

2zo-a c ' 

T^max — C^min)^ ^ 

2Z(j + CKn 


, OLr 


2zq + a c 

x [a 


1/2 


x \a 


min ? ^max ? 


mm? ^max 


(A.31) 


The conclusion follows as before using a max ^ 2zq and the symmetric relation —a m i n ^ 


2z 0 . 


□ 


Fix a triple (ao, ai, 02 ), with as usual the at increasing integers of sum N. We denote by 
Ao the common value of the a| — ay, where {/ky}ye{o, 1 . 2 } are given in (A.3). For arbitrary 
A E [—A c , A c ] we define the quantity 


£°(A) = ^ [a 0 a 0 (A) + aiay(A) + a 2 a 2 (A)] 


(A.32) 


where the ay (A) are three distinct roots of x 3 — x — A, numbered in such a way that 
ay(A 0 ) = ay. By construction, we have So(Ao) = 0 . 

Proposition A.l ensures the existence, for ( 7 , A) eD, of a continuous family x*(j, A) of 
solutions of (A.13), such that x*(0, A) has aj coordinates equal to ay(A). We set 


1 N 

(A) = ^£4(7, A). 


i= 1 


It follows directly from Lemma |A.2| that 

£o(A) - V7 ^ S 7 (A) ^ S 0 (A) + V7 • 


(A.33) 


(A.34) 


If A ey X 7 (A) changes sign in D at some A*( 7 ), then £*( 7 , A* ( 7 )) is indeed a stationary 
point of V/, satisfying the constraint 1 * £ S, Assuming for the moment that such a point 
exists, the following result characterises its signature. 


Lemma A.3. Assume that ( 7 , A) E IwtD', where D' C D is defined in (A.15). Then 
any stationary point x * of the family Bk, with triple (ao, ai, 02 ) = (0, M — k, M + k), is a 
local minimum of the constrained potential Vy. Furthermore, there exists a constant cq > 0 
such that if 7 ^ co\/A c — | A|, then any stationary point x * of the family Ck, with triple 
(ao, ai, 02 ) = (1, M — k — 1, M + k), is a saddle of index 1 of the constrained potential V/. 


Proof: First we note that by definition of D' , the function g defined in (A. 18) satisfies 


1 A 2.. .. 2 4 

“ W ~ f c ~ } ~ vs > vs s m “ ’ 


(A.35) 


which implies that points in V_ have a first coordinate x satisfying x < —l/\/3, and thus 
3x 2 > 1. By symmetry, points in V+ also have a first coordinate satisfying 3x 2 > 1. 
Stationary points x* in the family Bk have all coordinates in V±- since they are deforma¬ 
tions of points with all coordinates equal to a m i n or a: max . The Hessian matrix H 7 of the 
unconstrained potential at any stationary point x* defines the quadratic form 


N 

V l-y (V, HryV) = YX3(X*) 2 
i= 1 


N 


l ) y2 i + ~ Vi + i Y 


i=1 


(A.36) 
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This form is clearly positive definite if x* £ Bk, showing that x* is a local minimum of the 
unconstrained potential. Thus it is also a local minimum of the constrained potential. 

In the case where x* £ Cf, it has exactly one coordinate x in Vo, for which one easily 
checks that 3x 2 < 1. Thus Hq has exactly one negative eigenvalue, showing that for 7 = 0, x* 
1-saddle of the unconstrained potential. By Proposition |3.1[ x* is also a 1-saddle of the 


is a 


constrained potential, so that there exists a vector v £ S such that (v, Hqv) < 0. In fact, one 
can deduce from (A. 8 ) that the negative eigenvalue of Ho is bounded above by —c\ \J\ Q — | A| 
for a ci > 0 , while its other eigenvalues are bounded below by ci \J\ C — | A|. Since the second 
term in (A.36) has an ^-operator norm equal to 7 (it is a discrete Laplacian, diagonalisable 


by discrete Fourier transform), the Bauer-Fike theorem shows that x * remains a 1-saddle 
of the unconstrained potential as long as 7 < C 2 y / A c — | A| for some C 2 > 0. 

To show that this also holds for the constrained potential, we can use the fact that 
the eigenvectors of a perturbed matrix move by an amount controlled by the size of the 
perturbation (see for instance PI Thm. 4.1]). In this way, we obtain the existence of an 
orthogonal matrix 0 7 such that dOy = 0 7 — 1 has order 7 /y / A c — |A| and D 1 = OyHyOy 
is diagonal, with the same eigenvalues as Hy. It follows by Cauchy-Schwarz that 


(v, H. y v) = (OyV, DyOyV) 

= (v, DyV) + 2 (SOyV, DyV) + (SOyV, DySOyV) 




c 3\/ A c — | A | + 


C47 


y/X c ~ 


+ 


C57 


Ac - |A| 


(A.37) 


for constants 03 , 04,05 > 0. This shows that for y/(A c — |A|) sufficiently small, (v,Hyv) < 0 
and thus x* is a saddle of index at least 1 of the constrained system. However, the index 
cannot be larger than for the unconstrained system, so that if must equal 1 . □ 


Proof of Theorem 3.6. If we denote by ±A(y) the upper and lower boundaries of D , then a 
sufficient condition for S 7 to change sign is 


s o(A( 7 )) > y/l and H 0 (-A( 7 )) < -y/j 


(A.38) 


Without limiting the generality, we assume Ao ^ 0. Then the first of the two conditions is 
the more stringent one. For the family B *., using the fact that a m \ n 
near A c we obtain 




) - cQ + + 0(X C - A) (A.39) 


for some constant c > 0. Since we also have A(y) ^ A c (l — |y) = A c — y/S'j, inserting this 


in (A.38) yields the result. The case of the families Ck is similar, noting that the bound on 


7 in Lemma A.3 ensuring that they remain 1-saddles is fulfilled under the condition (3.6). 


In the case of the family Bo, one can obtain sharper bounds by first noting that x* has 
exactly half of its components in V_ and the other half in V+. Using the bounds given in 


Lemma A.2 we see that (A.34) can be strengthened to 


^o(A) — -\/7 ^ ^ ^o(A) + -\/7 


(A.40) 


Furthermore, we have 


A1q(A) — 2 


i(A) + -a E 


c(A) = ~2 «c(A) 


(A.41) 
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A sufficient condition for the stationary point to exist is thus 


-a c (A c (l-§ 7 )) > V7- ( a - 42 ) 

By definition of a c , this is equivalent to A c (l — § 7 ) > ^ 7(7 — 1). Taking the square yields 
the condition 27y 3 — 1357 2 + 54y — 4 < 0, which holds for 7 < | — \/5. □ 


B Proofs: Metastable hierarchy 

B.l Hierarchy of the B k 

Proof of Theorem \4-4\ When 7 = 0 , the value of the potential is constant on each family 


B k and C k . Using the expressions (A.3) of the atj, one obtains for these values 
aN(N - a) 


V 0 (B k ) = - 


4(N 2 - 3aN + 3a 2 ) ’ 


Vo{Ck+ 1 ) — - 


aN 2 — (a 2 + 8 a — 8 )N + 9 a(a — 1) 


4(1V 2 - 3aN + 3a 2 - 3a + 3) 


where a = M — k in both cases. Taking differences and simplifying yields 

(a — 1)(21V — 3a ) 3 


(B.l) 


V 0 (C k+1 ) - V 0 (B k ) = 
V 0 (C k ) - V 0 (B k ) = 


4(1V 2 - 3 aN + 3 a 2 ) (AT 2 - 3aN + 3a 2 - 3a + 3) 
(N - a - l)(3a - N) 3 


4 (N 2 - 3aA + 3a 2 )(N 2 - 3 (a + 1)JV + 3a 2 + 3a + 3) 


=: h\(a) , (B.2) 

=: / 12 (a) . 


Computing derivatives and proceeding in a similar way as in the proof of Theorem 3.5 


one 


obtains that a 1 —> hi (a) is decreasing, while a 1-7 /i 2 («) is increasing. Furthermore, it is 
immediate to check that h±(N/2 ) = h,2(N/2). We thus obtain the inequalities 

■ ■ ■ < Vo ( 6 * 2 ) —Vo (Bo) < Vo (Ci) —Vo (.Eh) < U 0 (di)-Uo(B 0 ) < U 0 (C 2 )-U 0 (Bi) < ... (B.3) 

(cf. Figurejb]). To prove ( |4.6| ), we have to check that relation (4.2) holds for each B k . Indeed, 
on the one hand we have 


(B h , U Bi) = V 0 (C k ) - V 0 (B k ) 
V i= 0 ' 


H[B k 


(B.4) 


for k ^ 2 , while on the other hand 
k 


h(b e , U Bi \ Bf] > Vb(Q) - Uo (Be) , i € [2, k - 1] , 

^ i =0 2 

(b q , U Bi) = Vo(Ci) - U 0 (Bo) • (B.5) 

V i= 1 / 


EE 


Thus the result follows from (B.3). 


When 7 > 0 is sufficiently small, the same partition still forms a metastable hierarchy, 
because the potential heights of the critical points depend continuously on 7 . □ 
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B.2 Hierarchy on B 0 

Proof of Proposition \f. 7[ The fact that allowed transitions between elements in Bq corre¬ 
spond to exchanging a particle and a hole follow directly from the connection rules (3.5) 


Indeed, any element in B\, with triple (0,M — 1,M + 1), is connected to M + 1 elements 
of i?o, which differ by a particle/hole transposition (see also Figure [4]). Any such transition 
affects at most 4 interfaces. Since the number of interface is always even, we obtain (4.12). 


In order to compute communication heights, we have to determine the heights of 1- 
saddles z* in C\. Recall that each of these saddles has 1 coordinate equal to a' 0l M — 1 
coordinates equal to a\ and M coordinates equal to a' 2 , where 

(a / 0 ,a / 1 ,a' 2 ) = ±a; / (l,l-M,M-2) , J = (M 2 — 3M + 3) -1/2 . (B.6) 


Plugging this into ( 4.10| ) yields 

F 7 (^(7))=Ro(z*(0)) 


+ l(u') 2 [M 2 4'/ai(^) + (M - 3 ) 2 I </a ' 2 (z*) + (2 M - 3 ) 2 I a ’ l/a ’(zn\ 


+ 0 ( 7 2 ), 


(B.7) 


where, similarly to (4.9), I a './ a 'X z *) denotes the number of interfaces of type a'/a'- of z*. 
The first-order correction to the height of the saddle thus only depends on the triple 

I(z*) = I a ' 0 / a i 2 {z*), , (B.8) 

where we use square brackets in order to avoid confusion with the triple (0 ,M — 1 , M + 1). 
Note in particular that I a ' 0 /a' 1 ( z *) + Ia' 0 /a' 2 ( z *) = 2, since only 1 component of z* is equal 
to «q. The following lemma allows to compare all these saddle heights. 

Lemma B.l. For all even p € [[2,1V — 2], the first-order terms V W in dBjj ) satisfy 

V^([2,0,p})>V^([0,2,p])>V^([l,l,p-l})>V^([2,0,p-2}) , (B.9) 

where [ a,b,c] stands for any saddle z* such that I(z*) = [a,b,c]. 


Proof: This follows from a straightforward computation, using (B.7) and the fact that 
(2M- 3)(M- 3) > 0. □ 

It remains to apply these expressions to the different transitions in Table [lj Consider 
for instance the transition shown in Figure [4j which is of type II.b. The two saddles 
encountered during the transition are of type [1,1 ,p — 1] and [0, 2,p], where p = 2 is the 
number of interfaces of the start configuration x*. Lemma [B. 11 shows that the second saddle 
is the highest. Combining this with the expression (4.10) of the height of x* yields 


H& = 


2(M — 3) 2 + p(2M — 3) 2 


4(M 2 - 3M + 3) ~ l> 4(M 2 - 3 M + 3) ’ 

which is precisely the expression given in the second line of Table [lj 


2(Af - 3) 2 - 3 p 


(B.10) 


The other cases are treated in a similar way. One just has to take care of the fact that 


the transition rules (3.5) allow for two possible paths, depending on whether ( 07 , 0 : 2 ) = 
(1,-1) + 0( 7 ) or (—1,1) + 0( 7 ). It is thus necessary to determine the minimum of the 
communication heights associated with these two paths. Table [2] shows the associated 
saddles in cases where the exchanged sites are not nearest neighbours. Table [3] shows the 
same when the exchanged sites are nearest neighbours. These saddle interface numbers are 
indeed those shown in Table DQ □ 
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1 

Il.a/b 

III 

IV.a/b 

V.a/b 

VI 

V 

( 0 , 0 ) 

(o,i)/(i,o) 

( 1 , 1 ) 

( 0 , 2 )/( 2 , 0 ) 

( 1 , 2 )/( 2 , 1 ) 

( 2 , 2 ) 

Hz*) 

[ 0 , 2 ,p + 2 ] 

[ 0,2 ,p] 

[ 1 , 1 ,P- 1 ] 

[ 0 , 2 ,p — 2 ] 

[ 0 , 2 ,p — 2 ] 

[ 0 , 2 ,p — 2 ] 


Table 2. Interface numbers of the highest saddle encountered along a minimal path between 
two configurations in Bq, if the two exchanged sites i and j are not nearest neighbours. The 
labels in the first row are the same as in Table [l] and rj denotes the number of interfaces of 
i and j. 



II.C 

IV.c/d 

V.c 

T] 

( 1 , 1 ) 

( 1 , 2 )/( 2 , 1 ) 

( 2 , 2 ) 

I(z *) 

[ 0,2 ,p\ 

[ 0 , 2 ,p — 2 ] 

[ 0 , 2 ,p — 2 ] 


Table 3. Interface numbers of the highest saddle encountered along a minimal path between 
two configurations in Uq, if the two exchanged sites i and j are nearest neighbours. 


Proof of Theorem 4-8. In a similar way as in the proof of Theorem |4.4[ , we prove that 
the relation (4.2) holds when the A p and A' p are ordered according to (4.15). Since all 
communication heights are the same when 7 = 0 , it will be sufficient to compare the first- 
order coefficients H^\ 

We start by showing that A 2 < A'±< ■■■ < A' M ,. For any p G [4, M' ], we note that 


H (1 \K. MuJ^U-U 4_ 2 ) = Ap) = 4(M^3M + 3 ^ • 


(B. 11) 


Indeed, the highest saddle encountered along a minimal path from A' p to A 2 U A' a U ... A ' p _ 2 
occurs during the type-III transition from A' to A p , and has interface number [ 1,1 ,p— 1 ], 
Since (B.ll) is a decreasing function of p, the condition (4.2) is indeed satisfied. 


Next we observe that 


H^{A 4 ,A 2 UA' 4 U---UA , M ,) = H^{A 4 ,A 2 ) = 


- 6 M 2 + 12 M - 12 


(B.12) 


4 (M 2 - 3 M + 3) ' 

Indeed, here the minimal path goes directly from A 4 to A 2 , via a saddle of type [0, 2, 2], The 
expression (B.12) is indeed smaller than (B.ll) for p = M ’, the numerator of the difference 
being bounded by AM 2 — 9 M + 12 which is always positive. 

Finally, we see that we have 

H {1) (A P1 A 2 U---UA' M ,UA 4 U---U A p - 2 ) = ff (1) ( Ap, a p - 2 ) = ~1^ 2 V^m + z) ’ (B - 13) 


the minimal path reaching communication height on a saddle of type [ 0 , 2 , p— 2 ]. Since (B. 13) 
is again a decreasing function of p, the claim follows. □ 


C Proofs: Spectral gap 

In order to prove Theorem |5.1[ we have to take into account the symmetries of the potential 
V~,. This will allow us to apply the theory in [5] on metastable processes that are invariant 
under a group of symmetries, which relies on Frobenius’ representation theory of finite 
groups (see for instance TO- 
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The potential is invariant under the three transformations 

r : (xi, .. .,x N ) ^ (x 2 , . ■ .,x N ,xi) , 

s : (xi,..., x N ) (x N ,..., xi) , (C.l) 

c : (xi,... ,xjv) ^ (~xi,-x N ) . 

It is thus invariant under the group G generated by these three transformations. This 
group can be written G = T>n x Z 2 , where £>at is the dihedral group of symmetries of 
a regular IV-gon generated by r and s, while Z 2 = {id, c} is the group generated by c, 
which commutes with r and s. The group G has order 4 N, and its elements can be written 
r l s^c k with i £ [[0, N — 1] and j, k £ {0,1}. It admits exactly 8 one-dimensional irreducible 
representations given by 

vrp<rr(rVc fe ) = p l a 3 r k , p, a, r = ±1 , (C.2) 

and N — 2 irreducible representation of dimension 2, whose characters are 

Xe,±(r l s j c k ) = TVt r A± (rVc fc ) = 2cos(^) 6 j 0 (±l) k , * € [1, f - 1] . (C.3) 

The basic idea of the approach given in [5j is that each of these N + 6 irreducible represen¬ 
tations provides an invariant subspace of the generator L of the Markovian jump process 
on Sq that approximates the dynamics of the diffusion. Thus the restriction of L to each of 
these subspaces yields a part of the spectrum of L. The eigenvalues of L can then be shown 
to be close to the exponentially small eigenvalues of C mm- We thus have to determine, 
for each irreducible representation, the smallest eigenvalue of —L. We do this in two main 
steps: first we show that the Arrhenius exponent of each smallest nonzero eigenvalue is 
given by the potential difference between certain stationary points in C\ and in A 2 , and 
then we compute the smallest prefactor of these eigenvalues. 


C.l Arrhenius exponent 


Each g £ G induces a permutation 7r 9 on the set of local minima So, leaving invariant each 
group orbit O a = {ga: g £ G}. Since V 1 is G-invariant, the generator L commutes with all 
these permutations. Thus there exist subspaces which are jointly invariant under L and all 
the 7 T g . Each of the irreducible representations of G provides one of these subspaces. 

Let 7r be one of the irreducible representations of G , and let d £ {1, 2} be its dimension. 
Then (5, Lemma 3.6] shows that the associated invariant subspace, when restricted to O a , 
has dimension da™, where 


a„ = 


I G a 


£ 

h^G a 


X(h) £ [0,i 


(C.l) 


Here x = T r7r denotes the character of n, and G a = {g £ G : ga = a} the stabiliser of a. We 
call active with respect to the irreducible representation 7r the orbits O a such that a™ > 0. 
Only active orbits will occur in the restriction of L to the invariant subspace associated 
with 7r; they are represented by a block of size da™ x da™. 

We select three representatives x* £ A 2 , y* £ B\ and z * £ Ci such that x* is connected 
to y* via z* in the transition graph Q. A possible choice is 


x (qu,...,qu,ou, a 2 ..., o?2) , 

z • • •, a 1, Uq, a 2 • • •, 0^2) , 

* / // // n // n\ 

y (^15 • • • ? ^15 ^2 5 ^2 • • • 5 ^2) 5 


31 




where a\ = 1 + 0 ( 7 ), 0:2 = —1 + 0 ( 7 ) and a ' 2 are each repeated M times, a\ and at" are 
repeated M — 1 times and a 7 } is repeated M + 1 times. The orbit 0 X of x * is precisely A 2 , 
and it has N elements. The orbits O y of y * and O z of z* have respectively 2N and 41V 
elements (they are proper subsets of B\ and C\). The associated stabilisers are given by 

G x = {id, r M s, r M c, sc} , 

G y = {id, r M ~ l s} , (C. 6 ) 

G 2 = {id} , 


where id denotes the identity of G and M = y. Note in particular that 

\G X \ . | G v 


I G x n G,, 


= 4 


I G x n Gy 


= 2 


(C.7) 


This means that each element in G x is connected to 4 elements in G y (via 4 saddles in G z ), 
and that each element in G y is connected to 2 elements in G x (cf. j5, (2.25)]). 

The possible Arrhenius exponents of eigenvalues of L are directly linked to which orbits 
are active for the different irreducible representations. We start with irreducible represen¬ 
tations of dimension 1, cf. (C.2). 


Proposition C.l. Let n be a 1-dimensional irreducible representation of G. Then 

• if M is even, then O x = A 2 is active if and only i/vr(s) = 7r(c) = 1; 

• if M is odd, then O x = A 2 is active if and only if tt(t) = 7r(s) = 7r(c); 

• if M is even, then O y is active if and only if n(r) = n(s); 

• if M is odd, then O y is active if and only if 7r(s) = 1. 


PROOF: The orbit O x is active if and only if ir(g) = 1 for all g € G x . Since 7r(id) = 1, 
ir(r M s) = 7t(?’) a/ 7t(s), ir(r M c) = 7r(r) A/ 7r(c) and 7r(sc) = 7r(s)7r(c), this holds if and only if 
7T(r) A/ = 7r(s) = vr(c). Similarly, the orbit O y is active if and only if it (?’) A:f_1 = vr(s). □ 


The corresponding result for the 2-dimensional irreducible representations given in (C.3) 
reads as follows. 


Proposition C.2. Let tt£± be a 2-dimensional irreducible representation. Then 

• O x = A 2 is active for 7T£ i+ if and only if l is even; 

• O x = A 2 is active for 7 r^_ if and only if l is odd; 

• O y is active for all representations ne t ±. 


Proof: By (C.3) we have x?,±(id) = 2, xe,±( rM s) = 0, Xt,±{ rMc ) = ±2cos(77r) and 
Xi,±(sc) = 0. Thus the sum (C.4) for O x is different from 0 for xi,+ if and only if t is even, 
and for xi,- if an d on ly if ^ is °dd. Since xi ,±(r M x s) = 0, the sum for O y is always equal 
to 1. □ 


Corollary C.3. The maximal Arrhenius exponent of all nonzero eigenvalues of the gener¬ 
ator is given by Vy(z*) — Vy(x*). 

PROOF: [5] Thm. 3.5] provides an algorithm determining the Arrhenius exponents for each 
irreducible representation 7r of dimension 1. They are obtained by replacing all inactive 
orbits by a cemetery state, which is at the bottom of the metastable hierarchy, and ordering 
all other orbits according to the usual hierarchy. If O x = A 2 is active and O y is inactive for 
7 r, then the largest communication height determining an Arrhenius exponent will be given 
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by V 1 {z*) — Vy(x*). If M is even, the representation it _|_+ has the required property, while 

if M is odd, this role is played by tt _ 

In the case of 2-dinrensional representations, 0 Thm. 3.9] shows that all communica¬ 
tion heights between active orbits yield Arrhenius exponents. The largest such exponent is 


obtained if O x and O y are both active, and Proposition C.2 shows that there are represen¬ 
tations for which this is the case. □ 


C.2 Eyring Kramers prefactor 

It remains to find the smallest prefactor associated with a transition of communication 
height Ty(z*) — Vy(x*). In the case of one-dimensional representations, [SJ Prop. 3.4] shows 


that the usual Eyring-Kramers law given in Theorem 4.3 has to be corrected by a factor 
\G X \/\G X n G y \ = 4 (cf. <[CT|). 


In the case of two-dimensional irreducible representations 7r^±, the relevant matrix ele¬ 
ments are given in [5] Prop. 3.7]. Alternatively, one can compute these elements “by hand” 
in the following way. We start by ordering the elements of the two orbits O x and O y 
according to 


O x = {x*, rx* 
Oy = { y*,ry * 


jv-i * 


X*}, 


,r^-y,r M q/V M+ V, 


The restriction of L to O x U O y consists in the four blocks 


Lxx — ^Qx^N i 

II 

-‘2q y i-2N , 



( 1 i 

(0) 

1 1 

(0)\ 

Lxy — Qx 






(0) 

1 

(o) •• 

• 1 


V 1 

1 

1 



,r 


t 1 
1 


M—l 


cy*} 


(C.8) 


Ly x — Qy 


( 0 ) 


( 0 ) 

1 

1 


\( 0 ) 


( 0 ) 


1 


l\ 


V 


(C.9) 


where 11,, denotes the n X n identity matrix, (0) stands for repeated zero entries, and 
|A_(z*)| /det V 2 E 7 i 


Qx — 


2 t r 


det V 2 E 7 (z* 


. e -[^(**)-^(**)]/e[i + 0 ( £ l/2| log£ |3/2jj ^ ( C . 10 ) 


while q y is a positive constant of order _ Equation (3.14) in [5] provides 

a set of vectors spanning the invariant subspaces associated with a given irreducible repre¬ 
sentation. Among these, we have to choose two linearly independent vectors for each orbit. 
A possible choice is 

u x = (2, x(r),..., xV^ 1 ), 0,..., 0) , 

u rx = {x{r N ~ l ), 2, xM> • • •, x(r N - 2 ), 0,..., 0) , 
u y = (0,..., 0, 2, x(r), • • •, X^ 1 ), 2, xM> ■ ■ ■, X^ -1 )) , 

u ry = (0,..., 0, x{ r ’ N ~ 1 )j 2, x(0> ■ ■ ■ > X^- 2 )^^” 1 ), 2, x(r),x(r N ~ 2 )) , (C.ll) 
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where x = X£,± is given by (C.3). In this basis, L takes the block form 


^xx ~ Qx 


4 0 
0 4 


L 


n _ n (2(x( r ) + 1) 

xy " x \ _o 


L yy ~ Qy 
L'ux Qy 


2 0 
0 2 

1 -1 

1 (x(r) + 1) 


(C.12) 


v -2 2 , 

We can now apply (5j Thm. 3.9], which states that the eigenvalues are equal to those of 

Ll x - Ll y (L; y )- l Ll x = -4 s in 2 (^jq x t 2 . (C.13) 


The result (5.1) follows, since the minimal value of the eigenvalues is reached for 1=1, 
both orbits are active for the representation 7Ti and this value is smaller than for all 
one-dimensional representations. 

Remark C.4. It is of course possible to obtain the same result directly from the expres¬ 
sions (3.15) and (3.17) given in [5] for the inner products (u,Lv), where u and v are basis 
vectors among (C.ll), even though these vectors are not orthogonal. It suffices to use the 
fact that the matrix elements of L can be obtained by computing 

-l 


/ 1 (u,v) \ 

(u,u) 

{u,v) 1 

\(v,v) / 


/(u,Lu) (u,Lv)\ 

(u , u) (u, u) 

{v , Lu) (v , Lv) 

\ (v,v) (v,v) ) 


(C.14) 


where for instance (u x ,u rx ) = cos(hr/M)(u x ,u x ). 0 


The last element of the proof of Theorem |5.1| is the following result on the Hessian 
matrices of Vo- 

Proposition C.5. The Hessian matrices of Vo at x * and z* satisfy 


det V 2 Vb(.x*) = 2 
detV 2 H 0 (z*) = - 


N -1 


M 


M —2 / 


(M — 3) M (2M — 3) 


2M—2 


(M 2 - 3 M + 3) 2M ~ 2 


= -2 


N—2 


[l + 0(N~ 1 )] , 


(M — 3)(2M — 3) 


A -' 2 ‘» = -lirat=- 1 + 0 T). 


(C. 15) 


Proof: We have already obtained invariant subspaces of the Hessian matrices in the proof 
of Proposition |3.1| However, since we used a non-isometric parametrisation of S, we cannot 
use expressions such as ( |A.7 ) directly to determine the eigenvalues. 

In the case of x* = (1,..., 1, — 1,..., — 1), it is sufficient to note that for any vector u of 
unit length in S, one has 


^2 Vo (x* + tu) 


M N 

= J> 2 C/"(1)+ u 2 U"(-1 ) = 2 , (C.16) 

0 i =1 i=M +1 


showing that in fact V 2 Vq(x*) = 211 ^-1, which has determinant 2 


JV-l 
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In the case of the saddle z * given by (C.5), the expressions (A.3) for the a'j yield 

U"{a\) = 3(a;) 2 - 1 = 


M(M - 3) 
M 2 - 3M + 3 
M (2M - 3) 
M 2 - 3M + 3 


= — 1 + 0{M~ 2 ) , 

= 2 + 3A7 -1 +0{M~ 2 ) 


U"(ol 2 ) = 3 ( c / 2 ) 2 - 1 = 3)(2 ^/ —^ = 2 - 3A/" 1 + 0(M“ 2 ) 

v 2; v 2 ' M 2 -3M + 3 v ' 


(C.17) 


We know from the proof of Proposition |3.1| that the (M — 2)-dimensional subspace of S given 
by Si = {aq + • • • + xm-i = 0, xm = • • • = xjv = 0} is invariant by the Hessian. Proceeding 
as in (C.16) with a unit vector u £ S\ shows that U"(a \) is an eigenvalue of V 2 Vo(-z*) of 
multiplicity M — 2. In an analogous way, the (M — l)-dimensional invariant subspace of S 
given by S 2 = 1x1 = ■ ■ ■ = xm = 0, xm+i + ■ • • + xn = 0} carries the eigenvalue U"(a' 2 ) 
with a multiplicity M — 1. This leaves a two-dimensional invariant subspace, for which we 
may choose the orthonormal basis given by the vectors 


v = 


w = 


V2M 




y/M(M- 1) 

The Hessian at (0,0) of the map (t, s ) Vq(z* + tv + sw) is found to be the matrix 

1 / 4M 2 - 15 M + 15 -3y/2(M- 1 )(M - 2) 

2 (A/ 2 - 3Af + 3) V-3 v / 2(AI- 1)(M - 2) -2(M 2 - 6M + 6) 

which has eigenvalues 

(M — 3)(2M — 3) 


and 


= -l + O^" 1 ) . 


2(M 2 - 3M + 3) 

The result follows via a Taylor expansion of log(— det V 2 Fq(-2 ; *))- 


(C.18) 

(C.19) 

(C.20) 

□ 
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